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■ Abstract We have performed measurements and numerical simulations on a bun- 

dle of vortex lines which is expanding along a rotating column of initially vortex- 

e*^ ■ free -^He-B. Expanding vortices form a propagating front: Within the front the su- 

perfluid is involved in rotation and behind the front the twisted vortex state forms, 
, which eventually relaxes to the equilibrium vortex state. We have measured the 

' magnitude of the twist and its relaxation rate as function of temperature above 

O ! 0.3 Tc- We also demonstrate that the integrity of the propagating vortex front re- 

^ O ^ ■ suits from axial superfluid flow, induced by the twist. 

PACS numbers: 67.57.Fg, 47.32.y, 67.40.Vs 

> 

' 1 Introduction 

00 ' 

' Since the pioneering works by Feynman, and by Hall and Vinen a rotating su- 

, perfluid has been associated with an array of rectilinear vortex lines, stretched 

• along the axis of rotation. This lowest energy state is known as the equilibrium 

I vortex state. Different kinds of collective excitations of vortex arrays have been 

. discussed over the years 1 1], including analogs of Kelvin waves on bundles of 

' vortices jHlsl] and Tkachenko waves HIH]. Recently a new state of vortex matter, 

the twisted vortex state, was experimentally identified and theoretically described 
in the B phase of superfluid ^He In this state vortices have helical configura- 
. , tion circling around the axis of rotation. Remarkably, certain such configurations, 

' uniform in the direction of the rotation axis, do not decay, although the energy of 

this state is higher than that of the equilibrium vortex array. The reason is that the 
force, acting on vortices, is directed along the vortex cores everywhere. 
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Fig. 1 (Color online) Numerical simulations of vortex dynamics demonstrate the twisted vor- 
tex state during the expansion of the vortex bundle along the rotating column. The subsequent 
relaxation of the twisted state proceeds from the top and bottom sample boundaries. The vor- 
tex configuration and the z-dependencies of the superfluid velocity components at given radial 
distance are shown at three points in time. The sample radius is 1.5 mm, length 40 mm, angular 
velocity 1 rad/s and temperature 0.5 Tc. 



The twisted vortex state is created when a bundle of vortex lines expands along 
an initially vortex-free rotating superfluid column, Fig. [ij Those segments of the 
vortex lines, which terminate on the side wall of the sample cylinder, propagate 
towards the vortex-free part and simultaneously precess around the central axis 
under the action of the Magnus and mutual friction forces. Such two-component 
motion leads to a vortex bundle which is helically twisted. Expansion along the 
column becomes slower with decreasing temperature as mutual friction decreases. 
Thus the spiral, created by the motion of the vortex end on the side wall, becomes 
tighter. One may expect that the resulting twist, characterized by the wave vector Q 
of the vortex helix in the bundle, becomes stronger as the temperature is reduced. 
We have measured the magnitude of the twist in the range 0.3 - 0.8 Tc and observed 
this expected behavior only at T > 0.45 Tc. Below 0.45 Tc the magnitude of the 
twist decreases again. 

In the real sample, which is not infinitely long, the twist cannot be completely 
uniform: At the top and bottom ends of the sample the vortices are perpendicular 
to the wall and the twist disappears there. The twist in the bulk unwinds when the 
vortex ends slide over the end plates of the sample cylinder. The model in Ref. |@] 
predicts that the relaxation of the twist becomes faster with decreasing temperature 
and distance to the wall. We have experimentally confirmed both properties and 
established reasonable agreement with the model. 

We also examine in this report the role of the twist-induced superflow in the 
propagation of the vortex front, which separates the vortex-free superfluid from 
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Fig. 2 Radial dependencies of the superfluid velocity from the simulation snapshot in the center 
panel of Fig.[T] Circles are the simulation results &i z= 1.8 cm, averaged over the azimuthal 
angle. Lines are fit to the twist model Eq. {T} with Q = 0.776/?^ ' . At r > 0. 1 2 cm the simulation 
results display the equilibrium vortex-free region around the central vortex cluster. This feature 
is ignored in the analytic model. 



the twisted vortex bundle. At T > 0.45 T^- the thickness of the vortex front in- 
creases while it propagates. At lower temperatures the twist-induced axial super- 
flow pushes vortices at the rear of the front forward. Eventually they catch up with 
the vortices in the head of the front and the front propagates in a thin steady-state 
configuration. 



2 Numerical simulations 

The essential features of the vortex front and the twisted state can be displayed 
by means of numerical calculations of vortex dynamics in a rotating cylinder. The 
simulation technique accounts fully for inter-vortex interaction and for the effect 
of solid walls 10] ■ In the initial state at f = the equilibrium number of vortices 
is placed as quarter-loops between the bottom and the cylindrical walls. During 
the subsequent evolution one observes the formation of the vortex front and a 
twisted cluster behind it, Fig. [T] The profile of the azimuthal component of the 
superfluid velocity Vs<;i shows an almost linear transition from the non-rotating 
state Vs0 = to equilibrium rotation Vs<;i w ^2/' within the region of the vortex 
front. This shear flow within the front is created by vortices which terminate on the 
cylindrical wall perpendicular to the axis of rotation. At the temperature of 0.5 
the thickness of the vortex front grows with time. Simulations at 0.4 on the other 
hand demonstrate a thin time-invariant front 1 6] . We discuss this difference below. 

The appearance of the twist is reflected in the axial superflow at the velocity 
Vsj, which is along the vortex expansion direction close to the cylindrical boundary 
and in the opposite direction close to the axis. As shown in Fig. |2] at constant 
height z the /--dependencies of v^i^, and Vsj are reasonably well described by the 
model suggested in Ref. lH]: 

(i2 + evo)r vo-Q^r' ... 

1 + Q^r^ 1 + Q^r^ 



where vq = (12 /Q)[Q^R^ l\n{\ + Q^R^) — 1] and/? is the sample radius. The wave 
vector Q of the twist has its maximum value close to the rear end of the front and 
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Fig. 3 (Color online) Experimental setup. 



decreases to zero at the bottom and top ends of the sample. The axial superflow 
affects the NMR spectrum of ^He-B and this allows us to observe the twisted 
vortex state in the experiment. 



3 Experiment 

The experimental techniques are similar to those in the recent studies of non- 
equilibrium vortex dynamics in -^He |8]. The sample of ^He-B at 29 bar pressure 
is contained in a cylindrical smooth-walled container with dimensions shown in 
Fig. [3] The sample is split in two independent B-phase volumes by an A phase 
layer, stabilized with applied magnetic field. In each B-phase volume an indepen- 
dent NMR spectrometer is used to monitor the vortex configuration. Two arrange- 
ments for pick-up coils have been used. They are labelled throughout this report as 
"setup 1" and "setup 2". These arrangements differ in the placement of the pick-up 
coils with respect to the upper and lower ends of the container, in Larmor frequen- 
cies, in the design of pick-up coils, and in the field homogeneity. The last property 
is mostly determined by the field distortion from the superconducting wire in the 
pick-up coils. For solenoidal coils in setup 1 it is AH /H w 6 ■ 10^^, while the coils 
in setup 2 distort the field more and AH /H w 1 .7 • 10^^. 

The initial vortex-free state is prepared by thermal cycling of the sample to 
temperatures above 0.1 Tc, where one waits at 12 = for the annihilation of the 
dynamic remanent vortices left over from the previous experiment [9]. After that 
the sample is cooled in rotation at 12 = 0.8rad/s to the target temperature. 

For vortex injection the angular velocity is ramped with the rate Q =0.03 rad/s^ 
above the critical velocity of the AB interface instability QcAB |1Q1 . after which 
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Fig. 4 (Color online) Examples of NMR records from a vortex injection experiment and defi- 
nition of the parameters Tcf, ^Lar and /itw The signal traces are recorded with the bottom spec- 
trometer in setup 2 at T = 0.50 Tc. 



Q. is kept constant, Fig. |4l In the instability event about 10 vortices are injected 
into the B-phase close to the AB interface. AiT < 0.6 Tc a turbulent burst immedi- 
ately follows and generates almost the equilibrium number of vortices 1 1 1]. These 
vortices then propagate towards the pick-up coil. The instability velocity X2cab is 
controlled by temperature and magnetic field profile and varies between 1 . 1 and 
1.5rad/s. A modification of this injection technique has been used at T > Q.^T^. 
Here the sample initially rotates at f2 > 2rad/s without the A phase. Then the mag- 
netic field is increased until the formation of the A phase starts. The AB interface 
immediately goes unstable and hundreds of vortices are generated [12]. 

While vortices expand along the sample column and pass through the pick- 
up coil their configuration is read using NMR. The signal is measured at a fixed 
frequency either at the location of the counterflow or the Larmor peak in the NMR 
spectrum jTll] . Absorption at the counterflow peak is sensitive to the difference 
in azimuthal flow velocities of the normal and superfluid components Vn^ — Vs^. 
It is at maximum in the initial vortex-free state with v^if, = 0. When increases 
the absorption decreases rapidly and drops to zero when v^ip ^ (l/2)vn0 fl3ll . 
Absorption in the Larmor peak is mostly sensitive to the axial flow v^z- It is zero 
in the vortex-free state while for a vortex cluster it has some finite value, which 
increases monotonously with increasing Vsz lfl3ll . 

An example of NMR measurement is presented in Fig. |4l To record signals at 
both Larmor and counterflow peaks with the same spectrometer the experiment 
is repeated twice in identical conditions. A rapid drop in the counterflow signal 
is seen which corresponds to the passage of the vortex front through the coil. We 
characterize this drop with time Tcf- The Larmor peak grows first to a maximum 
value htw, which corresponds to the maximum twist. Then the signal relaxes ex- 
ponentially with the time constant TLar to a value characteristic for an equilibrium 
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Fig. 5 (Color online) Strength of the twist as a function of temperature. The measurements are 
performed using the bottom spectrometer in both setups of Fig. [3] The ratio of the maximum 
amplitude of the Larmor peak in the twisted state to the amplitude of the Larmor peak in the 
nonrotating sample is plotted on the left axis. The maximum value of the twist wave vector Q, 
obtained in simulations, is plotted on the right axis. The solid curve shows the fit QR = 0.7(1 — 
a') I a. The broken curve shows the minimum magnitude of the twist for which propagation of 
the vortex front in a thin steady-state configuration is possible. 

vortex cluster. The propagation velocity Vf of the vortices at the head of the front 
can be determined from their flight time between the injection moment, where 
12 (f) = i2cAB. and the arrival to the edge of the peak-up coil, as monitored with 
the counterflow peak. We do not discuss Vf in detail in this report. 



4 Results 

The temperature dependence of the magnitude of the twist is presented in Fig. |5] 
From the experiment the raw data are plotted: The maximum height of the Larmor 
peak /ztw (Fig. |4]l, normalized to the height of the Larmor peak at 12 = at the 
same temperature. The data obtained in setup 1 and setup 2 are not identical since 
the spectrum shape around the Larmor peak depends on the magnitude of the 
magnetic field via the order parameter texture and also the homogeneity of the 
magnetic field is important: In setup 2 the homogeneity is worse and the Larmor 
peak cannot grow as high (i.e. as narrow) as in setup 1. Anyway, a clear maximum 
in the twist-induced signal is observed at around 0.45rc. The same behavior of the 
magnitude of the twist is confirmed in the numerical simulations: The maximum 
value of the twist wave vector behind the front, as determined from a fit of the 
velocity profiles to Eq. (|T]), also peaks at 0.45 Tg. 

The initial growth of the twist with decreasing temperature is expected. From 
the equations of motion of a single vortex which ends on the cylindrical wall, the 
estimate for the expansion velocity along z is 



Vu = a(r)[Vn0(/?)-Vs0 (/?)], 



(2) 
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where a is a mutual friction coefficient fl^. As the normal fluid velocity at 
the side wall Vn(i>{R) = ^2/? and the superfluid velocity induced by a single vor- 
tex can be neglected in our conditions, Eq. (|2]l gives vl; w a£2R. For the az- 
imuthal velocity of the vortex end in the rotating frame a similar estimation gives 
VL<j> « (1 — a')Q.R, where a' is another mutual friction coefficient lfl4ll . Thus 
the trajectory of the vortex end on the side wall is a spiral with wave vector 
Q = VL<j) /{Rvlz) « (1 — a')/{Ra). This value can be used as a rough estimation 
of the Q vector in the twisted vortex state. As one sees from Fig.|5] the simulation 
results at r > 0.45 indeed follow this dependence. 

For the decreasing twist at T < 0.45 Tc a couple of reasons can be suggested. 
First, the twist can relax through reconnections between vortices in the bundle, 
which become more prominent with decreasing temperature. Second, relaxation of 
the twist proceeds in diffusive manner within the twisted cluster. The source of the 
twist is at the vortex front, while the sink is at the end plate of the cylinder, where 
the twist vanishes because of the boundary conditions. The effective diffusion 
coefficient increases as the temperature decreases |@]. The faster diffusion limits 
the maximum twist in the finite-size sample at low temperatures. 

The relaxation of the twisted state is observed as a decay in the amplitude of 
the Larmor peak. The measured time constant TLar of the exponential decay is 
presented in Fig.|6]for three different positions of the detector coil with respect to 
the end plate of the sample cylinder. It is clear that the relaxation indeed becomes 
faster with decreasing temperature. It is also evident that the relaxation of the 
twist proceeds faster as the observation point moves closer to the end plate of the 
sample cylinder. Both of these features are well described by the model presented 
in Ref. [fi]. According to this model the relaxation time of the twist is 



Yi\r) [(l-a')2 + a2]' ^ ' 
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where C ^ 1 and b is the distance along z to the end plate, where the twist van- 
ishes. The fit to this model using C as the only fitting parameter demonstrates 
reasonable agreement with the measurements in Fig. |6] The experimental points 
show somewhat faster temperature dependence than the model. The discrepancy 
may result from shortcomings of both the model and the experiment. The model 
is constructed for the case of weak twist (QR ^ 1). In the experiment the probe is 
not point-like, but has a height comparable to the value of b itself. 

When Eq. (|2]l is applied to vortices within the vortex front, the following prob- 
lem arises: At the head of the front the superfluid component is almost at rest, 
Vsi|) ^ Vn0, and the expansion velocity of vortices is Vf w a£2R. Behind the front 
the density of vortices is close to the equilibrium and v^ij) w Vn,;, . Thus vortices at 
the tail of the front are barely able to expand and the thickness of the front should 
rapidly increase in time. On the other hand simulations show that vortices at the 
tail of the front do expand. Fig. [T] Moreover, at sufficiently low temperatures the 
expansion velocity of these vortices reaches the velocity of the foremost vortices, 
so that the front propagates in a thin steady-state configuration [6]. 

The explanation is that the vortex state behind the front is not the equilib- 
rium state, but the twisted vortex state. Taking into account the axial superflow, 
induced by the twist, Eq. (|2]l should be modified as vl^ = (X [v„(i,{R) — Vs^(7?)] + 
(1 — a')vs7(^). Given that v^ziR) is in the direction of the front propagation and 
Vs(|) (R) < Va(f, (/?) in the twisted state, the expansion velocity Vt of the vortices in 
the tail of the front is enhanced. This velocity can be estimated taking Vs-{R) and 
Vs(p{R) from Eq. ([T|: 



Vt = aQR 



1 



I -a' 1 



a QR 



1 ^ 



[l + Q^R^)log{l + Q^R^) 



(4) 



This velocity has a maximum as a function of the Q vector. If a/(l — a') > 0.52 
(i.e. T > 0.46 Tc lfl4h l the maximum value of Vt is less than the velocity of the 
foremost vortices Vf w aQR. In these conditions the thickness of the front in- 
creases while it propagates. When T < 0.46 Tc a wide range of Q values exists for 
which formally Vt ^ Vf. The minimum possible value of Q is shown in Fig.|5]by 
the broken curve. In these conditions the vortex front propagates in a steady-state 
"thin" configuration. These simple considerations become inapplicable, however, 
when vortices interact strongly within the thin layer. To determine, say, the sta- 
ble front thickness and the magnitude of the twist in the regime below 0.46 a 
different approach would be needed. 

The change in the front propagation at T « 0.45 Tc is observed in both exper- 
iment and simulations. In the experiment the decay time of the counterflow peak 
TcF in Fig. |4] can be used to extract the front thickness. The decay of the coun- 
terflow peak starts when the head of the vortex front arrives at that edge of the 
detector coil, which is closer to the injection point. The counterflow signal van- 
ishes when the last part of the front, which still possesses enough counterflow to 
generate the NMR response, leaves the far edge of the detector coil. The product 
TcpVf has the dimension of length and can be called the apparent thickness of the 
front. When the real thickness of the front grows with time its apparent thickness 
depends on the distance of the observation point from the injection point and on 
the rate with which the real thickness increases. When the front remains thin the 
value of its apparent thickness equals the height of the pick-up coil c/coii = 9 mm. 
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Fig. 7 Apparent thickness of the voitex front as function of temperature. The solid line is the 
prediction of the model in Eq. {5)- 



The measurements of the apparent thickness of the front are presented in 
Fig. |7] The scatter at higher temperatures has two sources. Partially it is due to 
the uncertainty in the determination of Vf from flight times which are not long 
compared to the measuring resolution. Another contribution is the variation in 
the decay profile of the counterflow signal which may display weak oscillations 
around a roughly linear decrease. At T > 0.45rc we have TcpV'f > c/coii and the 
apparent thickness decreases with temperature. Finally at T < 0.45 Tc the front 
becomes thin compared to the height of the pick-up coil. Assuming that at the 
moment of injection the front is infinitely thin we can write 



(5) 



where d is the distance from the injection point (i.e. position of the AB interface) 
to the nearest edge of the pick-up coil and V* is the expansion velocity at the 
position in the front where the NMR signal from the counterflow vanishes. Given 
that the latter condition roughly corresponds to Vs0 ^ ( 1 /2)vn0 we take V(* = (Vt + 
Vf) /2 if Vt < Vf and simply Vj* = Vf otherwise. Using Vt from Eq. Q and the simple 
estimates QR = (1 — a')/o: and Vf = aQ.R, we get from Eq. (|5]l the solid line in 
Fig.|7] which is in reasonable agreement with the experiment. 

In the simulations the thickness of the front grows with time at higher tempera- 
tures. This process slows down as the twist increases with decreasing temperature 
(Fig. |5]l. Finally at T w 0.45 Tc the twist reaches the value which is enough, ac- 
cording to Eq. to support a thin front configuration. At lower temperatures the 
thickness of the front indeed becomes time-independent and roughly equal to the 
radius of the sample. Simultaneously the twist behind the front starts to drop, as 
has been discussed above. It remains, however, within the limits where the twist- 
induced superflow is able to keep the front thin. 
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5 Conclusions 

We have studied the formation and relaxation of the twisted vortex state in ■'He- 
B in the temperature range between 0.3 Tc and 0.8 Tc. At higher temperatures 
T > 0.45 Tc the twist behind the propagating vortex front grows with decreasing 
temperature as Q (1 — a')/a. Here the thickness of the front increases while it 
propagates along the rotating column. The axial superflow, induced by the twisted 
state, boosts the expansion velocity of vortices in the tail of the front. This en- 
hancement increases as the twist grows with decreasing temperature. Finally at 
T w 0.45 Tc vortices in the tail of the front are able to catch up with vortices in the 
head. At lower temperatures the front propagates in a thin steady-state configura- 
tion, while the twist starts to decrease. 

The relaxation of the twist in a sample of finite length is not related to the 
front. It proceeds from the walls which limit the length of the sample along the 
rotation axis. The relaxation speeds up with decreasing temperature, unlike many 
other processes in vortex dynamics. 

While the understanding of the front propagation and of the formation of the 
twisted vortex state in the high- temperature regime T > 0.45 is quite good, in 
the low-temperature regime of a thin front the theoretical understanding is lacking. 
Especially interesting would be to consider the role of vortex recoimections and 
turbulence, both of which should increase with decreasing temperature. 
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